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Nonzero chemical potential studies with Wilson fermions should avoid the proliferation of flavor-equivalent 
nucleon states encountered with staggered formulation of fermions. However, conventional wisdom has been that 
finite baryon density calculations with Wilson fermions will be prohibitively expensive. We demonstrate that the 
volume method applied to Wilson fermions gives surprisingly stable results on a small number of configurations. It 
is pointed out that this method may be applied to any local or nonlocal gauge invariant quantity. Some illustrative 
results for < iptp > and < J > at various values of /t in a quenched lattice simulation are given. 



1. INTRODUCTION 

The behavior of quenched lattice QCD at fi- 
nite chemical potential is poorly understood. In 
previous simulations using staggered fermions an 
early onset of the baryon phase transition oc- 
curs at ^ opposed to the expected onset 
at where is the mass of the light- 
est baryon[l]. One possible explanation for the 
anomalous behavior of chemical potential is based 
on the stronger interaction than in nature of the 
4 degenerate quark fiavors present in the usual 
staggered formulation [2]. If this proposal is cor- 
rect, then going to a single light fiavor, as is the 
case for Wilson fermions, should substantially im- 
prove the situation. 

Previously, calculations of thermal expectation 
values (THEVs) in the Wilson case have been 
hampered by the large amount of computer time 
necessary to invert the quark matrix starting 
at all inital space-time points. However, recent 
investigations of disconnected quark loops have 
been made possible by global algorithms which 
either use stochastic estimation techniques [3] or 
exploit gauge in variance [4] to estimate fermion 
matrix elements. We investigate here the use of 
the method in Ref. 4 (which we refer to as the 
"volume method") in calculating THEVs at finite 
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chemical potential. We will see that this method 
gives numerically stable results on a surprisingly 
small number of gauge configurations for both lo- 
cal (< ipip >) and nonlocal (< J >) gauge invari- 
ant quantities. 

2. VOLUME METHOD 

The Wilson fermion action with chemical po- 
tential /J, is given by 

Sf= J2 (1) 

where the collective indices {i} and {j} include 
space-time, color and Dirac indices and where 

fc = l,2,3 

+ (1 + 7k)Ul{x - ak)6^^y+a^] 

-K[e~''{l - j4)U4{x)6^^y-a^ 

+e''(l + ji)ul(x-ai)6,^y+aA- (2) 

In a standard manner, one can calculate the 
THEV: 

J2 < ^(x,t)^(x,t) >= 

Y,{Tr[M-H£,t;£,mu. (3) 



0.0 0.1 0.2 0.3 



.5 0.6 0.7 



0.0 0.1 0.2 0.3 



.5 0.6 0.7 



Figure 1. Values of < tptp > as a function of jj, 
at K = .148. The horizontal line simply indicates 
the central value of < tptp > at = 0. The two 
sets of vertical bars represent the error bars on 
measured values of ^ and 



(The notation {•••}[/ on the rhs indicates a 
gauge field average.) Because of the average over 
gauge configurations in Eq.(3), we may replace 
the point-to-point loop propagator M~^(x,t;x,t) 
with a version which is summed over all inital 
space-time positions. That is 

< Mx,mx,t) >= Y,{TriM-H£,mu (4) 
where 

M-i(£",t) = ^M-i(£,t;£',t'). (5) 

This is because 

{Tr[M-\x,t;x',t')]}u = 

{Tr[M-\x,t;x,t)]}u6s,s'6t.t' (6) 

as a consequence of Elitzur's theorem[5]. The 
right hand side of Eq.(4) is simple to evaluate 
because it is just the trace of a quark propagator 
starting at all inital space-time points with unit 
weight. 



Figure 2. Values of the conserved charge density 
< J > as a function of at k = .148. The mean- 
ings of the vertical lines are the same as in Fig. 1. 



Likewise, the volume method can be applied to 
the THEV of any gauge invariant quantity. The 
conserved and appropriately normalized nonlocal 
Wilson charge density operator at nonzero /j, is 
given by 

J{x,t) = K{e^tl>{x,t + 1)(1 + j4)Ul{x,t)xp{x,t) 
-e-''xPix,t)il-ji)U4ix,t)xPix,t + l)}. (7) 

(The factors of e^^ will remove any trivial n de- 
pendence in < J > from quarks which do not loop 
in the time direction.) We then have the THEV: 

J2 < ■Hx,t) >= 

nJ^iTri-e" M-\x,t){l + j4)Ul{x,t) 

+e-''M-\x,t + 1)(1 - ji)lMx,t)]}u, (8) 

where we have again replaced the point-to- 
point space-time propagators by initial space- 
time summed ones, using Eq.(5). Since these ex- 
pressions assume a gauge field average, an im- 
portant numerical question is how well the gauge 
invariant signal in Eqs.(4) and (8) is projected 
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out of the gauge variant noise. For this purpose 
we turn to the numerical simulations. 

3. RESULTS 

We present results on a small number of config- 
urations in Figs. 1 and 2 for < tptp > and < J >, 
respectively. 

We used 5 gauge field configurations at each 
jj, (.3, .34, .4, .5, and .6) except at = and .36 
where 10 configurations were employed. Periodic 
boundary conditions were imposed on the quarks 
in the space direction and antiperiodic bound- 
ary conditions in the time direction, appropri- 
ate for finite temperature. The results come from 
16^ X 24 lattices at /3 = 6.0 and a relatively heavy 
quark mass hopping parameter of k = .148. (kc is 
approximately .1568.) In our conjugate gradient 
algorithm we found that it was very important to 
impose the convergence criterion directly on the 
quantity of interest. In this case we required that 
the relative change in < tptp > over the previous 
20 iterations to be less than 2 x 10~®. (The num- 
ber of iterations increased from about 60 at = 
to about 400 at = .6.) Our raw, unsubtracted 
result for < > at = was 11.2352 ± .0015. 
This is in agreement with 11. 23619± .00031 found 
with the Z2 stochastic method[6]. 

Fig. 1 represents the results of the < tptp > 
measurements as a function of jj, using Eq.(4). 
(Results in Figs. 1 and 2 have been normalized 
by the total number of space-time points.) The 
vertical lines on this figure (as well as the next) 
represent the error bar ranges on ^ and 
across 24 gauge field configurations (m^ = .689 ± 
.006, niA = 1.200 ± .022) [6] of which ours are a 
subset. (The A is the lowest mass single fiavor 
three quark state.) As can be seen, there is a 
statistically significant departure from the zero 
chemical potential THEV when /j, is below 

Fig. 2 shows the measurement of < J > as 
a function of /j, using Eq.(8). Again, there is a 
signicant departure from the vacuum value when 
/J, is below ^f^. (At = we measure < J >= 
.0004 ±^0015.) For iJ> the values begin to 
fall (< ipip >) or rise (< J >) sharply. 

As is evident from the figures, the volume 
method gives statistically stable results for both 



< ipip > and the charge density, < J >. This 
is very encouraging considering the small num- 
ber of configurations used. Thus, quenched finite 
density Wilson simulations are feasible without 
excessive amounts of computer time. 

Since there is only one light Wilson fermion fia- 
vor, our results are compatible either with a shift 
to lower /J, due to finite temperature [2] or a true 
early onset at = ^[7]. In order to clearly dis- 
tinguish between these two interpretations, fur- 
ther measurements at lower quark masses, where 
the pion and delta mass scales are more widely 
separated, are planned. 
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